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ABSTRACT 

The increasing weight and speed of carrier based aircraft are 
taxing the limit of conventional piston type energy absorbers which 
are used for arresting gear. Turbo-type absorbers have been proposed 
as an alternative. This study investigates a turbo-type energy 
absorber for an aircraft arresting gear under development by the 
Naval Air Engineering Center, A mathematical analysis, computer 
simulation and performance prediction are given for each mode of 
operation. It was initially expected that the absorber could operate 
in just two modes, forward and reverse, but the analysis shows that 
four distinct modes are possible, since not only shaft rotation but 
the path of fluid flow can be reversed. Theoretical performance 
predictions are also compared with test data from an existing smaller 
scale version of the absorber. The agreement is excellent. An 
approximation of scaling effects on power absorption is also 
presented. It is concluded that the turbo-type absorber is basically 
adequate to meet the demands of carrier operation for the forseeable 
future. 
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TABLE OE SYMBOLS 



Symbol 


FORTRAN 


Definition 


Units 


Latin 








a 


- 


Opening at blade exit 


ft 


b 


B 


Blade height 


ft 


c f 


CF 


Local skin friction coefficient 


— 




CL 


Pressure loss coefficient for 
lower passage 


- - - - 


c 

P 


CP 


Specific heat of fluid 


BTU/lb/°R 


C 

u 


CU 


Pressure loss coefficient for 
upper passage 


- - - - 


D h 





Hydraulic diameter 


ft 


h 





Static enthalpy 


ft-lb/ slug 


H 


H 


Total head 


ft-lb/slug 


B T 


- 


Total enthalpy 


ft-lb/slug 


1 R 


IR 


Rotor incidence angle 


degree 


X S 


IS 


Stator incidence angle 


degree 


S 


KB 


Blockage factor 


- 


l a 


LL 


Length of meridional contour of 
lower passage 


ft 


L 

XL 


LU 


Length of meridional contour of 
upper passage 


ft 


m 

m 


— 


Mass flow rate 


slugs/sec 


M 

f 


- - - - 


Frictional moment 


lb-ft 


P 


SP 


Static pressure 


lb/ft 2 


P E 





Equivalent relative total enthalpy 


O p 

ft / sec 
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Symbol 


FORTRAN 


Definition 


Units 


G> 


- - - - 


Power 


lb-ft/sec 


iSvef 


POWREF 


Reference power 


lb-ft/sec 


P T 


TP 


Total pressure 


lb/ft 2 


R 


R 


Radius 


ft 


r . 
1 


RI 


Radius of curvature of inner 
contour of passage 


in 


^Ne 





Effective Reynolds number 


— 




REL 


Reynolds number lower passage 


— 




RER 


Reynolds number rotor 


— 




RES 


Reynolds number stator 





^u 


REU 


Reynolds number upper passage 





XU 





Blade spacing 


ft 


r 

o 


— 


Radius of curvature of outer 
contour of passage 


in 


T 

R 


TORQUE 


Torque exerted on fluid by rotor 


lb-ft 


T 

Rref 


TORREF 


Reference torque 


lb-ft 


T 

S 





Torque exerted on fluid by stator 


lb-ft 


u 


— 


Internal energy 


ft-lb/slug 


U 





Peripheral velocity of blade 


ft/sec 


V 


V 


Absolute velocity 


ft/sec 


f 

V 


— 


Useful absolute velocity 


ft/sec 


V M 


VM 


Meridional component of 
absolute velocity 


ft/sec 


V 

u 


,vu 


Peripheral component of 
absolute velocity 


ft/sec 


W 


' w 


Relative velocity 


ft/sec 


t 

w 





Useful relative velocity 


ft/sec 
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Symbol 


FORTRAN 


Definition 


Units 


W. 

l 


- - - - 


Relative velocity lost due to 
incidence 


ft/sec 


¥ 

u 


VU 


Peripheral component of 
relative velocity 


ft/sec 


y r 


YR 


Rotor blade loss factor 





Y s 


YS 


Stator blade loss factor 





Greek 


O' 


ALPHA 


Absolute flow angle 


degree 




A3B 


Tangent to stator blade mean 
camber line at station 3 


degree 


a r 


AAR 


Aspect ratio of rotor 


----- 


A s 


ARS 


Aspect ratio of stator 





O' 

4B 


A4B 


Tangent to stator blade mean 
camber line at station 4 


r> r> r\ 


P 


BETA 


Relative flow angle 


degree 


P 1B 


BIB 


Tangent to rotor blade mean 
camber line at station 1 


degree 


P 2B 


B2B 


Tangent to rotor blade mean 






camber line at station 2 


degree 


A 


DEL 


Change 







ZETAL 


Channel bending loss coefficient 
for lower passage 


- 


^u 


ZETAU 


Channel bending loss ’ coefficient 
for upper passage 


- - - - 


^u 


ETAMU 


Friction efficiency for upper 






passage 


_ — « 


\e 


ETAML 


Friction efficiency for lower 






passage . 


•“ •“ •— *“ 




- 


Rotor efficiency 


- 
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Symbol 



X 

X 

X 



f 

l 



u 



"h 

Viff 

\ 



n 



0 ) 



p 

r T 



CD 

CD 



ref 



FORTRAN 

LMBD 

LAMDA1 

LAMDAR 

NR 

PIPOWR 

PITORQ 

PIOM- 

RHO 

TARR 

PHI 

PHI 

PSI 

OH 

OMREF 



Definition 

Wall friction factor in straight 
ducts 

Friction coefficient for lower 
passage 

Friction coefficient for upper 
passage 

Kinematic viscosity 
Dimensionless head 
Reference power 
Dimensionless Torque 
Reference rotational speed 
Fluid density 

Dimensionless torque coefficient 

Wall shear stress 

Flow coefficient used in Program 
NORMAL 

Convenient flow coefficient used 
in program REVERSE 

Velocity coefficient 

Rotational speed 

Reference rotational speed 



Rnits 



ft^/sec 

Ib-ft/sec 



slugs/ft'’ 
lb/ ft^ 



rad/sec 

rad/sec 
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Subscripts 



R Rotor 

S Stator 

Z Lower passage 

u Upper passage 

0 Outer passage 

1 Inner passage 

TH Theoretical 

MAX Maximum value 

ref Reference value 

1 Rotor inlet 

2 Rotor discharge 

3 Stator inlet 

4 Stator discharge 



Superscript 

t 

Useful component 
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I. INTRODUCTION 



Over the years, the speed, size, and weight of carrier based air- 
craft have steadily increased, thus demanding an ever larger energy 
absorption capacity from the shipboard arresting gear. Up to now, 
design improvements to meet these increasing demands have been limited 
to mere enlargement of the basic linear hydraulic energy absorber, 
which has been of the piston type. Modern aircraft, however, have 
loaded this system to its upper limit, and make any further development 
along these lines impractical due to weight and space limitations. The 
need exists for a new concept which will minimize space and weight 
demands while providing the necessary increase in energy absorption 
capacity. 

The solution proposed involves the use of a turbo-type rotary 
hydraulic energy absorber. In essence the fluid is put into motion by 
a suitably designed rotor and is retarded by a stator and by the 
resistance of the flow passage. The power input to the rotor is 
absorbed by the hydraulic resistance of the system, and is eventually 
dissipated as heat. 

The energy absorber is directly coupled to a cable storage drum. 

The drum profile is programmed to match drum rotational velocity to that 
required for optimum retarder performance. A device of this type was 
introduced in 15*68, and is known as the M-21 energy absorber. The 
initial design gives satisfactory performance, but is of rather limited 
capacity, thirty million ft. lbs. Therefore, a program has been under- 
taken to redesign the absorber to increase the capacity to fifty-five 
million ft. lbs. 
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A development program was initiated by the Naval Air Engineering 
Center (NAEC) (Ref, l). A full scale test of the programmed drum 
coupled with the M-21 absorber was conducted to determine the 
feasibility of the overall concept and to obtain performance data 
needed in the redesign program. A schematic diagram of the proposed 
system is shown in Figure 1 , 

A theoretical approximation of the flow conditions in a machine 
of this type was derived by Professor Michael H. Vavra (Ref, 2). To 
facilitate investigation of the effects of various design parameters 
NAEC formulated a computer program. From the computer analysis, a 
workable combination of design parameters was chosen such that torque 
and energy requirements were satisfied at the design point, the con- 
dition at which the absorber is expected to perform most often. 

Professor Vavra has also been commissioned to develop and test a 
flow model of the new design. This model will be useful for refining 
losses, checking for flow separation, and obtaining information that 
is not available from the mean streamline analysis. 

In this paper, a theoretical simulation and performance analysis 
of the energy absorber in all modes of operation is presented in terms 
of dimensionless parameters that facilitate comparisons of various 
designs and interpretation of model tests. The retraction problem 
is also studied, and preliminary estimates of the power needed to 
retract the programmed drum are given, 

A throttling system is incorporated into the design to aid in 
arresting aircraft of various sizes. The discussion and analysis 
presented here are concerned with the condition of open throttle only. 
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FIGURE I. PROGRAMMED DRUM ARRESTING GEAR 



II. DESIGN VALUES AND CONSIDERATIONS 



The main features of the energy absorber design that is analyzed 
here are summarized in Figures 2 and 3, and Table I. 

The customary approach in the design of turbomachinery is first to 
optimize conditions with respect to a single specified design point. 

The resulting performance over a broad range of off-design conditions 
is then investigated to determine overall acceptability. The condition 
at which the absorber is expected to operate most often is naturally 
chosen as the design point. Specifically, a design torque of 
340,000 ft. lbs. is required at a rotational speed of 76 rad. /sec. 

The maximum rotational speed the cable will impart to the drum, and 
consequently to the rotor, is about 85 rad. /'sec. This limit is fixed 
by the allowable landing speed of the aircraft, and by drum design. 

The mechanical power imparted to the rotor is continuously 
converted to thermal energy of the fluid. The temperature rise of the 
fluid is significant and imposes a constraint on the design. The 
temperature limit is the maximum allowable temperature the fluid may 
attain without causing cavitation. Obviously, heat must be removed to 
maintain equilibrium conditions. For this reason, provisions are made 
to bleed off a portion of the circulating fluid continuously and to 
replace it by fresh cool fluid. During a rapid series of arrests, the 
cooling problem is especially critical. 

The fluid used in the absorber is an ethyl glycol mixture. This 
fluid is acceptable because it remains usable over a wide range of 
temperatures. 
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SCHEMATIC DIAGRAM OF ENERGY ABSORBER 



FIGURE 2 



1 ft 





Figure 3 DIAGRAM OF ADSORBER BLADING 
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Table I PARAMETERS OP THE ABSORBER 



R.j = 13*^1 inches 
R^ = 16.25 inches 
R^ = 16.07 inches 
= 14.93 inches 
b,, = 2.02 inches 

It 



1! 

CO 

r£> 


2.02 


r = 

0 


4.18 


r . = 
1 


2.02 


(L . + 

' U1 


L 

uo 


2 




( L Ai + 


L Ao 


2 




Number 


of 1 


P 1B = 


57° 


P 2B = 


53° 


a 3B = 


57° 


a 4B = 


53° 


«B1 " 


• 95 


^B2 ~ 


.87 


^3 = 


.90 


^4 = 


.87 



= 10.26 inches 



= 10.26 inches 



72 



U). 



MAX = 85 rad, /sec. 

Fluid properties: 
p = 2.095 slugs/ft. 2 
v = 5*76 x 10”^ ft. 2 /sec 



c^ = 



c = 



.016 

.7 
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A detailed discussion of the effects of various parameters and loss 

coefficients is undertaken in the analysis section. A few words, 

however, are in order here. The radius of curvature r. of the inner 

1 

contour of both the upper and lower passage must not be smaller than 
the blade height. This insures that the flow does not separate, and 
it also increases absorption capacity. Estimates of losses in this 
machine are obtained through various empirical turbomachinery cor- 
relations as is shown in the analysis. Enthalpy losses in the passages 
are significant. Much of this loss is due to the bending of the 
passage with friction responsible for the remainder. However, a 
change in the coefficient of friction c^ which expresses the direct 
wall friction losses in the passages has little influence on overall 
performance. A value of c^ = 0.016 was chosen since it gives best 
agreement with M-21 test data. 

The blade shapes first optimized by Professor Vavra from a 
performance standpoint were found by NAEC to be unacceptable on a 
stress basis. Blade thickness, spacing and angles were accordingly 
revised by NAEC, and a blade of sufficient strength was found; its 
characteristics are listed in Table 1, 

The reliability of this analysis can only be judged experimentally 
by correlation with test data. The theory correlates very well with 
actual test results on the M-21 design. From this it appears that the 
theoretical analysis is reasonable and that the new design should be 
workable, if not optimum. 

Naval carrier operations often require a series of rapid and 
repeated arrests. The drum must be retracted thus reversing the 
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rotational direction of the rotor. This process is investigated to 
obtain preliminary estimates of the power required to retract the 
absorber. 
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III. ANALYSIS OP THE PROBLEM 



A. MODES OP OPERATION 

The energy absorber is capable of operation in four modes. These 
modes are summarized in Table II. 



Table II MODES OP OPERATION 



Plow Direction 


Shaft 


Rotation 




Normal 


Reverse 


Normal 


Mode 1 


Mode 2 


Reverse 


Mode 3 


Mode 4 



Mode 1 corresponds to design point conditions. In this mode, the 
fluid passes successively through stations 1, 2, 3» 4 and back to 1, 
in that order. Note that the flow leaving the rotor at station 2 is 
moving radially outward. See Figures 2 and 3. For fluid circulation 
in this sense, conditions at stations 2 and 4 are particularly sig- 
nificant. At the rotor exit, station 2, the direction of the flow 
with respect to the rotor, is assumed to remain constant at all flow 
rates. At stator exit, station 4 > the direction of the flow with 
respect to the stator is likewise assumed to remain constant. While 
this assumption of fixed relative and absolute flow directions is not 
exact, it represents an excellent approximation when the blades are 
closely spaced as in the present design. 

In Mode 2, the sense of the fluid circulation is exactly as 
described above. However, the direction of shaft rotation is reversed 
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Naturally, stations 2 and 4 are still governing in the sense that flow 
directions at these stations are known. 

In Mode 3 > the shaft rotation is the same as for Mode 1, hut the 
direction of fluid circulation is reversed. Consequently, the fluid 
passes successively through stations 2, 1, 4* 3 and back to 1, in 
that order. Note that the rotor exit is now at station 1 , and the 
direction of flow at this station is radially inward. In this case, 
stations 1 and 3 > are controlling. The relative flow direction at 
rotor exit, station 1, and the absolute flow direction at stator exit, 
station 3» are now assumed to remain fixed. 

In Mode 4» "the sense of fluid circulation is the same as in 
Mode 3 i tut the direction of shaft rotation is reversed. Of course, 
stations 1 and 3 still govern flow direction as in the previous case. 

The same set of equations describe Modes 'I and 2, but a separate 
set of equations is needed to describe the reverse flow Modes, 3 and 4 
Note that only two sets of equations are needed for the four Modes, 
and each set is governed by the flow direction relative to the blade 
exits. Two separate computer programs are needed since one program 
will be used to solve each set of equations. Program Normal finds the 
solution for Modes 1 and 2, and program Reverse finds the solution for 
Modes 3 and 4 * Both programs will be discussed in the appendix. All 
quantities and results discussed throughout the paper are obtained 
from these programs. 

The velocity diagrams for each mode are presented in Figures 4 » 5 » 
6, and 7, All angles and velocities are drawn to scale. 

The thermodynamic processes in the absorber are shown on the 
enthalpy-entropy diagrams, Figures 8 and 9« The process for normal 
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Scale; 0. 5 inch - 1 00 ft/sec 
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Figure 4 Velocity Diagrams for Mode 1 Operation at 

„ “max R 2 , 

= " 4 ' 8 X 10 
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Figure 5 Velocity Diagrams for Mode 2 Operation at 



^Ne 



“MAX R 2 
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4.8 x 10 
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Figure 6 Velocity Diagrams for Mode 3 Operation at 

^e “ " 4*3 * 10 6 
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Figure 7 Velocity Diagrams for Mode 4 Operation at 
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Figure 8 Total Enthalpy-Entropy Diagram of the Thermodynamic Process 
in the Absorber for Normal Circulation, Modes 1 and 2 
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Figure 9 Total Enthalpy-Entropy Diagraja of the Thermodynamic Process 
in the Absorber for Reverse Circulation. Modes 3 and A 



circulation, Modes 1 or 2, is shown in Figure 8, while that for 
reverse circulation, Modes 3 or 4> is shown in Figure 9* The 
direction of shaft rotation affects the exact quantitative enthalpy 
values. The graphs, however, are merely qualitative . 

A "total head" H (per unit mass) is used which differs from the 
standard total enthalpy H^,. This is discussed in greater depth in the 
analysis section, but a few words are needed at this point to clarify 
the discussion. 

The fluid is considered incompressible which allows the thermal 
and mechanical effects to be uncoupled. A "total head" is 
conveniently defined as 




Henceforth H is referred to simply as head. The head and total 
enthalpy are linked by the simple relation 

H = - u 

where u represents internal energy per unit mass. 

B. DIMENSIONAL ANALYSIS 

Results are presented in terms, of appropriate dimensionless 
quantities. The use of non-dimensional terms simplifies the cor- 
relation of data, the comparison of various designs, and the inter- 
pretation of model tests. 

The three quantities chosen as consistent dimensional reference 
parameters are the radius of the rotor outlet R£, the fluid density P, 
and the shaft rotational speed ou or ou^. Y , whichever is appropriate. 
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Using these parameters, one defines the effective Reynolds number 
for this design to be 

wR 2 2 

*N " — 
e 

where v is the kinematic viscosity. And the dimensionless total head is 




H 

w 2 R 2 2 



In the discussion of results, it is convenient to define the flow 
coefficient $ in the following specific way, namely, 




where the algebraic signs of cu and $ are always positive for normal 
shaft rotation and negative for reversed shaft rotation. This sign 
convention for shaft rotation is maintained regardless of the direction 
of the fluid circulation, hence the use of the absolute magnitude | | 
in the definition. 

In the actual computer programs NORMAL and REVERSE, however, some- 
what different sign conventions are more convenient. All equations are 
written in such a way that the three quantities V^ 2 wR^ and § are always 
positive. For this reason the four modes are computed as separate cases 
For comparison of various modes and designs the dimensionless torque 

power and rotational speed coefficients are respectively 
T 

® MAX R 2 P • 
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n 6> = 



<? 



“^MAX R 2 5p 



n 



U) 



gj 

0) 

MAX 



The detailed analysis considered later also uses various 
dimensionless quantities for convenience in the derivations. The 
parameters used to facilitate the non-dimens ionalizat ion in that case 
are not necessarily the same as introduced in the present section, and 
the separate cases should not he confused. 



C. ENERGY BALANCE. AT EQUILIBRIUM 

In general the rotor adds energy to the fluid. In traveling 
through the passages and returning to its initial state the fluid 
dissipates this mechanical energy through various losses. Under 
equilibrium these losses exactly cancel the original energy input. 

It turns out that the energy input of the rotor can be expressed 
in terms of the flow coefficient § in the form 

A % = A + B * 

where A and B are true constants. The losses can be expressed as 

““loss - A ' + B ’ * + c ’ ** 

» , t 

where A , B , and C may be considered constant to a first approximation. 
Actually these quantities are weak functions of $, a fact which must 
be considered in the exact solution. The energy relations which fix 
the equilibrium state are shown in Figures 10 and 11. Figure 10 
represents the two modes which correspond to normal flow direction. 
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of Normal Circulation 



Mode 4 I Mode 




J?3 



Figure 11 shows the two modes for reverse flow. These results are 
drawn to scale, except that the parabolic lines are only approximate. 

The final equilibrium points, where the losses balance the energy 
input, however, and the relative distribution of losses at these 
points are in correct scale. 

Initially, it was thought that only two of the above four modes 
would prove to be physically possible, one corresponding to normal 
shaft rotation, the other corresponding to retraction. Surprisingly 
enough, the present analysis has demonstrated that operation in all 
four of the modes is theoretically possible. This point is clearly 
indicated in Figures 10 and 11. 

D. MODEL TESTING 

Professor Vavra has proposed thax a model testing program be under- 
taken. The model will not be a true dynamic model, but will be a purely 
static device that uses air as the working fluid. Construction of the 
model will constrain the air flow to the normal direction only. However, 
the scale model should be useful in checking for flow separation 
through the passages and refining the loss coefficients' for the 
analysis of the absorber. The effects of the tangential component of 
velocity will also be investigated. 

The discussion of the above static model is not within the scope 
of this thesis. However, a hypothetical model is considered which is 
a true dynamic model of the actual full scale device. The size of the 
hypothetical model is taken as identical to that of the static model. 

The model is scaled to 0.569 the size of the actual absorber. In 
both cases the fluid is air. When expressed in dimensionless terms, 
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the performance of the hypothetical model is basically the same as that 

of the full scale machine, except for Reynolds number effects. The 

full scale machine operates at a maximum Reynolds number, as defined 

above, of 4*88 x 10^ while the maximum Reynolds number of the model is 

5 

3,70 x 10 . The results for this model are presented along with the 
results for the real absorber, 

E. ASSUMPTIONS AND LIMITATIONS OP THE ANALYSIS 

A one-dimensional mean streamline analysis is obviously only a 
crude approximation of the actual performance of any turbomachine, 
but is a starting point necessary in most development programs. The 
validity of the assumptions used are very important to the reliability 
of the analysis. The primary assumptions used in this x’eport are 
summarized below. 

The fluid is treated as incompressible. This assumption is 
excellent as long as the provisions to maintian the fluid at a 
temperature lower than the cavitation temperature are adequate. This 
conveniently uncouples mechanical and thermal effects and allows a 
form of "total head" H to be defined which is independent of 
internal energy. 

If the velocity distribution across every section were turly 

uniform, the meridional component V at. the mean streamline could be 

7PH 

computed simply from a continuity relation of the form 

where Q is the volumetric flow rate and A,, is the meridional cross- 
section area. The actual velocity distribution, however, is ncn- 
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uniform because of the effect of streamline curvature, boundary layers 
and so on. To allow for these effects, approximately, the actual 
meridional velocity is assumed to be governed by a relation of the 
form 

v m-*b4 

where, the so-called blockage factor Kg is treated as a constant. The 
blockage factors were extrapolated from M-21 test data. Their reliability 
in the reverse flow modes is doubtful, however, since all test data 
were obtained only from Mode 1 operation. An extension of these data 
to all modes would be very useful. 

The present analysis treats operation of the absorber as a steady 
state phenomenon. Strictly speaking, the operation is actually 
unsteady. However, the arrest and retraction processes require a 
sufficient number of drum revolutions that the steady state assumption 
provides a fair approximation. 

A major assumption is that the relative flow direction at the rotor 
outlet and the absolute flow direction at the stator outlet are 
constant. For the case of normal flow this is an excellent ap- 
proximation, However, the validity is reduced for the reverse flow 
modes. In Modes 5 and 4 > "the flow encounters a sharp leading edge, 
and exits past a blunt trailing edge. The implication is that the 
flow is not as well guided as in the normal flow case. Thus the 
actual turning angle of the flow is undoubtedly less than predicted 
by the theoretical analysis. The power absorbed in both reverse flow 
modes would consequently be less than predicted. 
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Losses, in general, are approximated by various empirical cor- 
relations and are reasonable as shown by test results. They are 
discussed in depth in the detailed analysis section. 

Incidence losses, however, pose a problem. The incidence loss 
approximations introduced in the analysis section are reasonable as 
long as the incidence angles remain small. This is obviously not 
the case for the retraction modes as shown in Figures 5 and 1 * The 
approximations break down physically at these severe incidence angles. 

It was thought earlier that two of the four theoretical modes 
might sho w themselves to be physically impossible. For example, a 
large incidence loss might conceivably preclude the attainment of the 
necessary energy balance. This has not proved to be the case, however. 

It appears from the present analysis that operation in all four modes 
remains possible physically. It is true, of course, that the predicted 
performance for the two retraction modes, Modes 2 and 4> will be 
inaccurate because of these large incidence effects, but this in itself 
does not eliminate the existance of these modes. 

At extreme angles of incidence, a distinct possibility of stalling 
exists. In normal operation, this limits the attainable power absorption. 
In a retraction mode, however, stalling would not be totally undesirable. 
If the rotor stalls, the energy absorbed by the device is greatly 
reduced, and the power necessary to retract the rotor is minimal. 

All of the foregoing assumtions are thought to be reasonable. 
Refinements of various points would be useful, especially for large 
incidence angles and blunt trailing edges. 
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IV. RESULTS AMD DISCUSSION 



The discussion of the problem up to this point has already included 
a number of significant results reached in the analysis as shown in 
Figures 4 through 1 1 . 

Additional detailed results for all modes of the absorber are 
presented in Figures 12, 13, 14, and in Table III. The results for 
Mode 1 of the dynamic model are also included. 

The main results of this investigation are summarized concisely in 
Figure 12. This is a logarithmic plot of the dimensionless power n 
as a function of the dimensionless shaft speed Performance 

of the proposed absorber, of the hypothetical model, and of the 
existing M-21 machine are compared on this graph. Note that the 
theoretically predicted performance of the M-21 . version as shown by 
the continuous line is in excellent agreement with the actual measured 
performance as shown by the data points. 

The actual horsepower is plotted as a function of absolute shaft 
speed for all four modes of the proposed absorber in Figure 13- The 
power absorbed during .arrest and the power expended during retraction 
are presented for a range of possible operating speeds. 

The maximum dimensionless power is also shown as a function of 

Reynolds number in Figure 14* This graph, although based on only 

e 

three computed points, does correctly portray, approximately, how 
scaling affects power absorption. It would be a simple matter to 
calculate additional points on this curve from program NORMAL. 
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POWER COEFFICIENT 




Figure 12 Hon— dimensional Power vs. lion— dimensional 
Rotational Speed 
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HORSEPOWER 




Figure 13 Absolute Power vs. Shaft Speed for the Four 

Operating Modes of the Proposed Energy Absorber 
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Table III. Continued 
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As mentioned earlier, the actual existence of four stable modes is 
a surprising result. The original expectation, based largely on 
intuition rather than analysis, was that only one mode would be possible 
for each direction of shaft rotation. But the absorber is an unusual 
machine which tends to confound intuitions based on experience with more 
conventional devices. Consider for example the fact that losses are 
necessary for the machine to absorb the mechanical energy input. The 
smaller the loss coefficients, the higher the final equilibrium flow 
rate will be, and the greater the energy absorption. Thus we are led 
to the seeming paradox that the lower the loss coefficients the higher 
the energy lost! 

The power curves of Figures 12 and 13 are interesting and 
significant in that all curves have nearly constant and nearly equal 
slopes. The power can be expressed in the form 



where the slope of the curve is roughly T]~3 for a ll cases. 

The horsepower curves of Figure 13 indicate that Mode 1, the 
design point mode, absorbs the greatest amount of energy. Mode 3> 
however, is a very close second. Refinement of the analysis should 
lower the predicted absorption of “Mode 3 somewhat. Nevertheless, the 
favorable comparison of this mode to the design mode is interesting. 
Figure 13 also shows that the power needed to retract the rotor is 
reasonably small for Modes 2 and 4» especially in comparison to the 
power absorbed during arrest. 

The question of which mode the absorber will actually operate in 
under various conditions of use is interesting, but by no means 
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completely clear at this time. It is thought possibly that if the 
absorber be started from rest with shaft rotation in the normal 
direction, Mode 1 should occur. If retraction follows immediately, 
the retraction mode will probably be Mode 2. This is based on the 
fact that during arrest the flow will build up a large momentum which 
will tend to persist in the same sense during the subsequent retraction. 
If retraction be started after the fluid motion ceased, however, the 
resulting mode of operation remains uncertain. 

Some method of sensing and/or controlling flow direction would be 
useful to insure that Mode t was occurring on all arrests. Mode 3 
operation might not have enough absorption capacity to handle large 
aircraft . 

The general question of which mode occurs is an area for further 
study. Testing of the actual absorber with the idea of definitely 
determining the operating modes under various conditions would be 
extremely useful. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



The development of the problem and the exact solution of the 
two sets of equations describing normal and reverse flow conditions 
clearly show that four stable operating modes exist. The analysis 
adequately predicts the absorption capacity of the proposed device 
for Mode 1. However, the analysis is not as accurate for the 
reverse flow or retraction modes. 

It is recommended that the retraction problem be further studied 
to refine loss coefficients and to more accurately predict, the power 
necessary to retract the absorber. It is also recommended that 
further study be made into the question of predicting which mode 
occurs under various operating conditions. This study would 
encompass the unsteady phase of the problem. Actual testing of the 
absorber should also be conducted with the intention of determining 
which mode will occur. A method of flow control or sensing would 
also be useful to guarantee proper performance in that Mode 3 
absorption capacity may be insufficient in certain cases. 

This turbo-type of energy absorber appears to be completely 
adequate to meet the demand of carrier operations now and in the 
forseeable future. The potential for future development and 
enlargement is ample. 
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APPENDIX A 



DETAILED ANALYSIS OF THE ENERGY ABSORBER WITH THE 
FLUID CIRCULATING IN THE NORMAL DIRECTION, 

MODES 1 AND 2. 

A one dimensional mean streamline analysis after Vavra (Ref. 2) 
will be developed with the following assumptions: 

- constant 
o ' = constant 

Uniform flow at stations 1 to 4 except for blockage factors. 

Blockage factors are independent of the flow. 

The fluid is incompressible. 

Steady state conditions exist. 

1.0 Torque and Velocity Relations 

The torque T_ exerted on the fluid by the rotor is positive 
K 

when acting in the direction of rotation. 

T R = “ [E 2 V u 2 - R 1 V . <U 

where m is the mass flow rate, 
m = p A V 

m 

= p (2tt R 1 b R ) V Mi (2) 

= p(2rr R 2 b R ) Kg^ 

The quantity 2tt R b R is the flow area, and V ^ - are the meridional 

velocities at the mean streamline, and and 

which are assumed to be independent of the flow. 



*12 



are blockage i actors 
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( 3 ) 



By introducing a flow coefficient 




all necessary velocity and thrust equations can he written in terms of 
known geometric parameters and the flow coefficient. An iteration 
technique is later employed in program NORMAL to calculate $, The 
following velocity relations are developed now in order to simplify 
the energy relationships which follow. 

The sign convention used is that all peripheral velocity com- 
ponents are positive when iij the direction of rotation. Prom Figures 
4 or 5 simple geometric relations are obtained which become, using 
eq. 3. 



W 2 = V M2 /cos 3 2 = cuR 2 $/cos P 2 



( 4 ) 



W u2 = V M2 Tan $ 2 = tuR 2 § Tan P 2 



( 5 ) 



V u2 = W u2 + tuR 2 = wR 2 ^ $ Tan P 2 ^ 



( 6 ) 



Also 



V u2 = V M2 Tan a 2 = U)R 2 § Tan a 2 



■ (?) 



An expression for or in terms of the known relative exit angle can be 
obtained by equating the previous two equations, or 



Tan o? 2 = Tan P 2 + l/$ 



( 8 ) 



similarly at the rotor inlet 




( 9 ) 




( 10 ) 
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The mass flow equation and the definition of the flow 



coefficient yield 



and 



R 2 

V M1 



S? 

Si 



S 

V = $ U)R — - 
ul 2 R 



S2 

- — Tan p + u® 

Si 1 1 



( 11 ) 



( 12 ) 



The mass flow equations must he extended to include flow 
through the stator. 



m - 2rr R 2 h R p V M2 Kg 2 2rr R^ h s p 

= 2n R 4 b s p V M4 S 4 



(13) 



From the above relations, the meridional velocity equations are developed. 



V 

Ml 


R 2 


S2 




V 

M2 


' R 1 


Si 




V 

M3 


R 2 


S2 


b R 


v 

142 


= S 


S3 


b 

s 


V 

M4 . 


R 2 


^2 


s 


V 

M2 


' R 4 




b 

s 



( 14 ) 



( 15 ) 

(16) 



1 , 1 Frictional Losses in Upper Passage 

Let be the frictional moment from forces along the walls 
of the upper passage from station 2 to 3, acting on the flow opposite 
to the direction of rotation, 

; H 3 V u 5 * “ E 2 \2 - M fu < 17) 
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Setting 



( 18 ) 



M- = \ m R V 0 
fu u . 2 u2 



and ^ = (1 - \) 
eq. 17 becomes 



R, 



V u3 _1 We, V u2 

3 



Prom the mass flow equation 




§ cuR 



R 2 ^~B2 

2R 3 *33 




b 

s 



09 ) 



( 20 ) 



( 21 ) 



1,2 Frictional Losses in Lower Passage 

With being the frictional moment caused by forces along 

the lower passage walls from station 4 to 1 acting on the flow opposite 
to the direction of rotation, the equations for the lower passage 
are developed in a manner similar to the previous development. 





= X. m R. V . 
A 4 u4 


(22) 




= 1 - 


h 




(23) 






R. 






V u1 


= \ie 


V — 

V u4 H 1 




(24) 


V M4 


= § CO 


H h. 

2 R 4 


^B2 b R 

V b s 


(25) 



From Figure 4 it can be seen that 



V . = V„. Tan ot 
u4 M4 4 



Using the previous three equations gives 

R 2 b R 

V u1 = % * " R 2 R^ ^ r Ta " *4 



( 26 ) 



(27) 
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From Figures 4 or 5» one can relate the absolute inlet angle 



to absolute velocity components as 

V 

Tan cv = ~~ (28) 

Ml 

Using the relations developed for the inlet velocity eqs. 11 and 27, 
the absolute inlet angle is 



R 



Si 



Tan a. 



'““i-’WFIb, *“ -4 



Similarly the relative inlet angle is 



Tan 8 1 = 



V „ - co R. 
ul 1 



Ml 



tuR . 



= Tan or - ^ 



Ml 



iitut ing 
fan 8 = 



yields 




Kbi m 

- Tan. 4 


- 4" fe 



(29) 



(30) 



The relative and absolute velocities can now be obtained from 
their respective components by employing the following general equations: 
„2 



2 2 
V M u 


(51) 


2 2 
= V M +W u 


(32) 


= V M /cos 8 


(33) 



52 



1 . 3 Torque Relations 

The torque as defined in eq, 1 can now be expressed in 
terms of known geometric parameters and the flow coefficient. Using 
eqs. 37, 6, and 12 in eq. 1 yields 



5 b R 2 

T R = 2TT V ^ P ^2 ® 






l Tsui fJ 2 - , W r vJ (34) 



The non-dimensional torque coefficient is now introduced as 

it 



T 



t R " 



R 



2n(R ) 5 ^ p K^ 2 co 2 
^ b 

s 



4 + }2 P 2 - V 



R 

~ Tan o' 
s 



4 



(35) 



The torque T g which the stator exerts on the fluid is positive 
when acting in the direction opposite to rotation 



T = m (R-; V , - R. ’ V .) 
s 3 u3 4 u4 / 



(36) 



Substituting the values for the whirl components gives 



T s - 2nR 2 5 Uf) p hi\u f (rj T£m “4 

Again introducing a non-dimensional torque coefficient 
T 



' b R \ ^B2 



T s = 



„ „ 5 ,' b R\ 2 _ _ 

2TT R 2 Wj P W ^B2 



(37) 



p • ^RP 

- V 4 + 4 LV Ts " p 2 - r Tsm “ 4 J 



53 



The general torque equation is 




The torque exerted on the rotor by the fluid minus frictional moment 
losses is equal to torque exerted on the fluid by the stator. 

The dimensionless torque coefficient introduced earlier is 
again given as 



2.0 Energy Relations 

The torque and velocity equations developed in the preceeding 
section remain in terms of the unknown flow coefficient §. In this 
section an algebraic energy equation will be formulated and solved 
to yield I. 

For the sake of this development, a total head term H will be 
defined which differs from the standard definition of total enthalpy 
H™ which is 




T 



(38) 



- h + ~ 



( 39 ) 



or 




(40) 



where h is the static enthalpy 



p is the static pressure 



u is the internal energy 



and is the total pressure. 
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Since the fluid is treated as incompressible, one can uncouple 
the thermal effects from the mechanical effects, and define a ’total 
head’ H which does not include thermal effects. This new term is 
denoted by 





( 41 ) 



Obviously the enthalpies are related by 



(Hij - u) = H (42) 

This new definition of total head H is useful in that constant H and 
constant lines are coincident. Throughout the remainder of the 
analysis, useful total head H will be referred to as simply head. Note 
that this total head has the units of enthalpy (per unit mass). 

The process undergone in the absorber is now represented in 
Figure 8, Obviously the process does not return to the original total 
enthalpy , although it does return to the original total pressure 
and head. 

From Figure 8 the energy rise across the rotor is shown to be 
equal to the energy losses across the stator and through the passages, 

or 

- AH u - AH^ - AH S = 0 (43) 

Therefore , 

(H 2 - - (H 2 - H 3 ) - (H 3 - H 4 ) - (H 4 - Hp = 0 (44) 

where all the above terms are written as positive quantities. 
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The four terms of eq. 44 will be derived in terms of the flow 
coefficient and the known geometric parameters as in Section 1.0. Eq. 
44 is then straightforwardly solved to yield §. [Program NORMAL is 
employed to accomplish this task. Since various loss parameters 
depend on an iteration technique must be used. 

2.1 Conditions in Rotor 

The incidence angle of the rotor blading i is defined as 

.K 



R 






IB 



( 45 ) 



or merely the angle between the tangent to the mean camber line at the 
leading edge and the relative inlet velocity. The incidence angle i 

K 

is seen to become more negative in the direction corresponding to a 
greater lift coefficient for the blade. 

The relative velocity can be resolved into two components 
respectively tangential and normal to the mean camber line. The energy 



associated with the tangential component W is assumed to be largely 
recoverable. The energy associated with the normal component is 
assumed to be totally lost. 

The useful relative velocity at the rotor inlet is then 



I 

= N, cos i 

± ± S\ 



( 46 ) 



and the velocity component lost due to incidence is similarly 



¥ u = w i sin S 



(•* 7 ) 



The basic energy equation of turbo-machinery is 
*12 • % = U 2 V u2 “ U 1 V ul 



( 48 ) 
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which upon expansion becomes 



^2 " ^1 “ 



2 2 2 
V +U -W 
2 2 2 



2 2 2 

v -hj -w 



(* 9 ) 



Using static enthalpies, regrouping, and cancelling terms in eq. 48 
gives 

2 2 2 2 

V - U 2 V - U 1 

h 2 + = h 2 + (50) 



Substituting for the static enthalpy and rearranging yields 

2 2 2 2 
p 2 \ P 1 W 1 U 2 " U 



2 — — " ( u 2 " 



( 51 ) 



Upon resolving W^/2 into its components and regrouping further, one 



obtains 



a + !A a . M i' g . v _l - U A / 

\ P 2 J \ p 2 2 / V U 2 " U 1 " 2 ) 



(52) 



This leads naturally to the definition of the "equivalent relative 
total pressure" as 

^ = <») 



2 ‘ 



and 



El 

P 



= 

V D 2 



2 2 

ir 2 - u . ^ 






( 54 ) 



and to the definition of the rotor blading losses as 

Ap f r, , w i in r 
— = |> 2 - up - —J = {- 



r> - P 
'El E2 



) 



( 55 ) 
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It also proves useful to express this loss in terms of a dimensionless 
coefficient defined as 



R 






El " P E2\ 



JE W " 
2 2 



(56) 



Recalling that 



H = H - u 
T 



(42) 



and using the basic energy equation, eq. 48, one obtains the energy 
rise across the rotor in the form 



(H - H. ) = U 0 V . - U. V . - {sin 2 i 

' 2 1 2 • u2 1 ul l 



2 2 

V V 

R 2 + Y R 2 



} 



(57) 



This equation can be conveniently non-dimensionalized by 
2 

dividing by U . The velocity expressions developed in Section 1.0 can 
not be utilized to simplify eq. 57, and the dimensionless energy rise 
across the rotor (II A H) can be written as 

H - H (sin 2 i ) . R .2 

(Wr ■ -hr 1 - 1 - — s- 2 - ( rt ) 

CU R 2 ^ 



{ { Tan P 2 - ' b7 1 W Tan ■%) COa2i R} 



(58) 



! 2 lIk^J + (j^- b s T w. Tan <0 ] 



- $ 



R 



2cos 2 |3, 



-} 
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where the dimensionless head coefficient is defined as 



n = -AiL_ 

AH 2 2 

u) R^ 



( 59 ) 



The loss coefficient Y of eqs. and 58 is a function of the 



deflection A3, where 



A 3 R = 3 2 - Pi 



and the aspect ration defined by 





(60) 



This will be further discussed in Section 3*0. 

2.2 Performance of Rotor Blading 

The rotor is designed for the maximum useful enthalpy rise 
across the blades. If this were a frictionless process, the useful 
energy rise would naturally be greater than in an actual process 
involving friction. Therefore, the rotor efficiency is defined as 



H 2 ~ H 1 

^ H 2 Th " H i 



(61) 



where the total heads H^, Hg, and are shown on the H-s diagram. 

Figure 8. Of course the pressure rise through the rotor is related 
very simply to the corresponding energy rise according to the 
relation 




P 



(62) 
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The numerator of eq. 6l is obtained from either eqs. 57 or 
58. The denominator of eq. 6l is similarly obtained, but for the 
theoretical case both the incidence and blading losses are non-existent. 
It follows from eq. 58 that 



^ ^Th^R 



r H 2Th ■ H 11 

12 2 J 

u) R 2 

1 + i {Tan p 2 - ^ 

s 



S4 



Tan 



} 



(63) 



is the theoretical dimensionless energy rise across the rotor. The 
efficiency can obviously be ■written as 



( n AH^r 

% “ <" "Wr 



(64) 



The general definition 



we: 



O. 

V l -1 o 



= m A H 



(65) 



It is convenient to introduce a dimensionless power coefficient as 




P 

3 v 5 
ptU MAX R 2 



(66) 



Recalling that 



V M2 
cu Rg 



m = 2u p R 2 b R 



*32 



(3) 

( 2 ) 



and 



n A H 



A H 

2_ 2 
tu R 2 



(59) 
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the dimensionless power and head coefficients can now be related as 



{ 2 "(£> *8 } 5 " A H 



( 67 ) 



MAX 



The power is related to the torque by the simple expression 



n ' ?= ^ y 



( 68 ) 



2.3 Energy Losses in the Upper Passage 

The energy loss in the upper passage, assumed as only a total 
pressure loss due to friction between the rotor and stator, is expressed 
as 



AH = H 0 - H = 

u 2 3 p 



T2 



T3 



(69) 



or 






(70) 



In a frictionless process 



H = ( h ) = !3_Ei + II 

2 K 3'T H 2 2 



(71) 



and from eq. 17 with the frictional moment M„ equal to zero 

p IU 

• R„ V _ 

(y ) = V — = — ^2- 

v u3'TH u2 R 2 

5 1 Mu 



(72) 



Resolving /2 into its components, one can write eq. 71 as 



p v V ‘ 

H ^ _ P 3TH . M3 , 1 U3 

\ ~ H 3 TH “ + "2“ + 2 



T1 2 Mu 



(73) 
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Similarly for an actual process 



H 



3 




( 74 ) 



the energy loss can now "be written using the two previous equations as 




(75) 



A convenient dimensionless pressure loss coefficient can now 
be denoted as 




(76) 



r 

where C will be shown to depend primarily on the radius ratio — of 
the upper passage in Figure 2. This is similar to energy loss 
coefficients developed for bends in smooth pipes. 

Eq. 75 is non-dimensionalized, and using eq. 76 one obtains 



(n ah) u 



c 

u 

2 



M3 



2 

U) ft- 




V 



u3 



2 

id IL 



(77) 



Using previously developed expressions for one obtains an equation 
similar in form to eq. 58 or 
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2.4 Performance of Stator Blading 



From Figure 4 the stator incidence angle i is 

s 



i = a - a? B 

s 3 3 



( 79 ) 



and the useful inlet velocity V 7 is 

3 



V, = V, cos i 

3 3s 



(80) 



The total head ahead of stator blades is given by 



3 p 2 p 



(81) 



while the useful head is 



H ' 

3 p 2 



T 3 



(82) 



The above three equations combine to yield 



, (V,') 2 - (V,) 2 

h 3 - H J + 2 2 



V, 2 

= - -j- (sin i o ) 



(83) 



Similar to rotor blading losses, stator blading losses will be defined 
as 



Y = 



P T3 ~ P T4 H 3 ~ H 4 



S P y 2 
2 V 4 



(V 4 2 /2) 



(84) 



or 



V 9 V 

= H 3 " Y s ~t = ( sin ‘ O Y 



iZ 

s 2 



(85) 
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Subtracting from both sides yields 

2 2 

H 4 - H 1 - H 3 - H 1 - < s “ 2i s> i" - Y s \ («> 

The same non-dimens ionali zing technique used in Sections 2.1, and 2.3 
is again employed and the energy loss in the lower passage is now 
available in the general form as 



H 4 - H 1 H 3 - H 1 4 , V, 2 

-» 2 r/ Vr/ ' 2 V 



(sin 2 i s ) 



+ $ [(sin 2 i s ) 2 Tan 3 2 j + ^ |(sirA g ) \_\jT~ 



2 i / . 2. 



1 R 



^2 



l-lu b s ^3 



j + Tan 2 3 2 






1 l\ 



2 l_R, b cos a.j j i 

h J 



v 4 s - 4 



(87) 



The above equation is actually 



h 4 ~ H 1 H 3 - H 1 H 3 ~ H 4 



2^ 2 
<u R 2 



2 d 2 
(0 R 2 



2 _ 2 
(0 R 2 



This can be further reduced to give 



h 4 - h 1 h 2 - Hl H 2 - h 3 h 5 - h 4 



2.2 
io R 2 



2 _ 2 
to R 2 



2 _ 2 
to R 2 



2 ' 2 
a) R 2 



( 88 ) 



which gives the energy losses in the lower passage in terms of the 
energy changes throughout the remainder of the system. The next 
step is to find the energy loss in the lower passage and solve eq. 

44. 



64 



B 








Again the loss coefficient Y is a function of deflection and 

s 

the aspect ration, where 



A(3 = a - a 

s $ 4 



(89) 



and 



A = 



s R, - R. 
3 4 



(90) 



2.5 Losses in Lower Passage 

The development of the energy loss in the lower passage 
closely follows the technique of Section 2.3. 

The energy loss in the lower passage is 



p 4 V 

“ H 4 " H 1 = f + 






(91) 



for a frictionless process 



R. 

(V ) = Y — =■ 

V Vl'Th u4 R 



(92) 



or relating theoretical to actual velocities yields 

V 



< V u1 ^Th 



ul 

\]A' 



(93) 



Therefore, 



2 

\ ^ P 1^Th V M1 , 1 , v u1\ 
H 4 - ^pTh ~ P + 2 + 2 It y/ 



while 



P-i V ^ v 2 

H„ = _L + IM__ + _ifl — 

p 2 2 



1 



(94) 



(95) 
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the pressure loss coefficient is 



( p 1^Th " P 1 



= C 



Ml 



i 2 



( 96 ) 



The expressions for and are now substituted into eq. 1 , 
along with eq. 6 to yield 

„ 2 



H 4 - H 1 " C f 



Ml 



-R., 



+ ^ ” 71 w) IrJ v - 



■u4 



(97) 



2 r 

This is easily put into the dimensionless form used previously or 

< naH b - t (r) + 0 - V 2 )(r r; Tan a S'\ (98) 



*B4 



Eq, 97 can now be substituted into eq, 88 and $ can be solved 
for straightforwardly. Iteration is necessary, however, since all loss 
coefficients depend on $, 

3,0 Loss Coefficient in Rotor 

Vavra ( 3 ) gives a relation for the losses in axial turbine 
bladings as 

-f “ °' 99 - A B “ ■ m 



where Ap is the flow deflection angle in degrees. The quantity Y is 
the velocity coefficient defined as the ratio of the actual and 
theoretical discharge velocities, which for the rotor is 



Y 



W„ 



R W. 



2Th 



( 100 ) 



and 

a P r = P 2 - P 1 ' (101) 
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The theoretical case occurs for Y^ equal zero which in 

It 

accordance to eq, 56 gives P^. = P„ orra , Substituting eq. 100 into 
eq. 53 yields 



P P T) W 

El _ E2Th _ ^2 1 _2 

p " P _ p + 2 ^2 



(102) 



The actual rotor blading losses using eq, 102 can now be written as 




(103) 



In accordance with Horlock (Ref. 4 ) "the loss coefficients depend 

ratio A- of eq. 60. Hence let Y., 

it a. 

of eq, 103 now denoted as Yp be the loss coefficient across the rotor 

S 

blading for A_ = 3 and R_ = lO”'. The latter term is deiinea as 
K IMit 



on Reynolds numbers R^ and aspect 






D hR W 2 



(104) 



where D, is the hydraulic diameter of the discharge area of the flow 

ila. 

channel between adjacent blades. If/414 is the rotor blade suacings at 
the radius and ap is the opening at the discharge as shown in 
Figure 3» an approximation of the discharge opening is given by 



Hl “Ai 003 P 2 

then 

_ 4 ( flow area) — 7 —^ 

hR - wetted perimeter - 2 (a^ + bp) 



(105) 



(106) 
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Substituting for and rearranging yields 



2b 



D. 



hR 



1 + 



_R 

\ 

^ cos S 2 



(107) 



The blade spacing .4^ can be found with respect to the radial 
extension R 2 - R^ by using the blade loading criterion of Zweifel (Ref. 5) * 
For axial bladings this approximation is 



,4U 



R 



0-425 



R 2 “ R 1 (Tan $ 2 - Tan P 1B ) cos 2 3 2 



(108) 



Substituting this equation into the expression for the hydraulic 
diameter yields 



D. 



hR 

R 



4 ) 



(b p /R 9 ) (Tan8 9 - Tanb^cosb 



(109) 



2 , + ' R' 2 



0,425(1-r 1 /r 2 ) 



Placing the above equation into the expression for the Reynolds number 
and recalling that 



W„ = 



iuR 2 $ 



2 cos 8, 



(4) 



one obtains 






§ tu R r 






cosg. 



(bfi/^) (TanPg - Tan8 1B ) cos 8 2 

0.425 (1 - R-j / r 2 



•] 



(110) 
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The loss coefficient is dependent on the aspect ratio and 
Reynolds number through the following relations after Horlock (Ref. 4 ). 



* 



R 




( 111 ) 




( 112 ) 



and 

y r - y r Gy (H3) 

t 

where the coefficient holds for the actual aspect ration only at a 

5 

Reynolds number equal to 10 . By combining the above equations the 
generalized loss coefficient is simply 




I?) - 1} (^y 25 



‘R 



t 

\ 



114) 



These relations allow the actual loss coefficient of the energy absorber 
to be determined under varied conditions# 

Obviously the Reynolds number is not totally specified by 

geometric parameters, but depends on the flow coefficient §# In 

1 

iterating for $, the process must be assumed to begin at Y , or a 

5 

Reynolds number of 10 , and subsequently correct the loss coefficient 
as improved values of § and R^ are calculated. This procedure must 
also be employed for the velocity coefficient Y which is dependent on 
$ through a dependence on the actual inlet angle 3^ . The iteration 
must be continued until convergence is achieved for all values 
dependent on $. 
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3*1 Loss Coefficient in Stator 

This section follows closely the development in Section 3*0* 
Eq. 99 is again used with 

V, 



Y = - 
s V 



4Th 



(115) 



and 



Ag = O' - a. 

s 3 4 



( 116 ) 



For a frictionless process in the stator, the total pressure 
("T4^Th s 1' a, 1' or ex i't would equal total useful inlet pressure P^ , 

or 



(Pm/) 



(V/)' 



T4 y Th _ T 3 _ i Th 



(117) 



Substituting eq. 115 into eq* 117 yields 

P p V 2 

111 m 14. + _1_ 4 

P P y ^ 2 

s 

For an actual flow 



Pm/ Pz V„ 
H _M = 14 + _4_ 

4 P P 2 



(118) 



(119) 



which yields through eq. 84 



1 * 
Y = - 1 = Y o 

s y 2 s 

s 



(120) 



This was as expected from Section 3*0* The loss coefficient can now be 
written after eqs. 112, 11 3 > and 1 1 4 as 



v ■ i (•>» * Y> - 1 



( 121 ) 
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i / 10 5y 2 5 
Y = Y (“■ — ) 

s 3 V V 



(122) 



and finally 



^ ♦'PM (£ 



5\* 2 5 



(123) 



where the aspect ratio is defined by eq. 90 as 

b 



A = 



s R-, - R. 

3 4 



The Reynolds number for the stator is 



D, Y, 



(124) 



where for a g “^coso^, the hydraulic diameter is 



2b 



D. = 
hs 



1 + 



b I & 
s' s 

COS O', 



(125) 



The blade spacing according to Zweifel (Ref. 5) is 



A' 



R 3 R 4 (| Tan - Tan Ijcos^o^ 



0.425 



( 126 ) 



and from eqs. 25 and 26 



R ~ b " ^2 



1 



V 4 = * “ R 2 R 4 \ V cos °4 



(127) 



The Reynolds number can now be expressed using the above four equations 



as 
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,b 

15 



iTl fS. 6 



5 iu R r 






V v 1 ~ 

b s 2 

cos + j- 1 - ( | Tan or^ - Tan o^-g | ) cos o/^ 



0.425 (R 5 /R 2 ) (l - 



( 128 ) 



A series of iterations is again needed to insure proper 
convergence of all parameters, 

4,0 Pressure Loss Coefficients in Upper and Lower Passages 

In accordance with eq. 76 the pressure losses in the upper 
passage are expressed by 



(Px)fPV, “ 



V, 



M3 



3 y Th 3 _ p 

p " u 2 

and similarly for the lower passage by eq. % as 



(129) 



kl>Th - P 1 



= C 



jM 1 

i 2 



(130) 



With Apb being the pressure drop in the bend, the loss 
coefficients of these bends defined as 



s = 



Apb 

fv 2 



( 151 ) 



are shown in Figure 15 as functions of the Reynolds number R^, where 
similar to eqs. 104 and 124 



V* D. 



(132) 



Vavra (Ref. 2) has developed the set of curves shown in Figure 15* 

These curves are based on the data of Sprenger (Ref. 6), Interpolation 
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was to obtain values for ratios of r^/b different than those found 
in Sprenger. 

The hydraulic diameter is again 



I) 



h 



4 flow area 

wetted perimeter 



( 133 ) 



For a ring channel with flow area of 2n E b the wetted perimeter is 
2(2ttR), or 




Hence at the rotor inlet there is 



\ = 2b H 

and at the stator inlet 



D h = 



2b 



(134) 



(135) 



(136) 



These above two values are introduced into eq. 132, with V 
equal to the corresponding value of V^, that is, before the rotor 

and Vj^ before the stator. The result is 



"Su 




(137) 






V M1 2b 



R 



( 138 ) 



The Reynolds number is then used to obtain the loss coefficients § of 
Figure 15* 

The loss coefficients § of Figure 15 have been obtained from the 
measured pressure drops by subratcting the pressure drops that ’would 
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occur in a straight duct at the corresponding Reynolds number and 
surface roughness. 

There are several conclusions which can be drawn from Figure 15 
and from additional data in Sprenger (Ref. 6). In general, the loss 
coefficients for bends made up of two concentric radii increase with 
decreasing r i/b ratios. The curves § versus R^ that exhibit a 
minimum indicate the flow separations occur at the inner radius, and 
should be avoided. Finally, Sprenger (Ref. 6) shows that two bends 
with r i/b = 0,5 arranged in series to form a 180° bend have a loss 

tz 

coefficient of 1.66^ at R^ = ICr where £ is the loss coefficient 
of one of the 99° bends. This data will be a first approximation in 
choosing the loss coefficient for the 180° bends in the energy- 
absorber. 



4 . 1 Estimate of the Pressure Loss Coefficients in the Upper 
and Lower Passage 

The losses in the upper and lower passages will be assumed 
due to two components, passage bending and friction. The losses due 
to duct bending were discussed in the previous section. An arbitrary 
assumption is to take these loss coefficients equal to 2.0 and 2.5 
for the upper and lower passage respectively. The values of £ are 
given in Figure 15 for various value of r i/b. It is recommended that 
r i/b be chosen between 0.5 and 1.0. Any bend with a value of r i/b 
larger than 1.0 will be assumed to have the same losses as a bend with 
r i/b equal 1.0. 

To the above losses must be added the frictional losses in a 
straight duct of hydraulic diameter and average length L, where 

L. + L 



L = 



(139) 
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and ]T is the distance along the inner contour and L q that along the 
outer contour in the meridional plane from the trailing edges of one 
blade row to the leading edges of the other row. (See Figure 2) 

In accordance with eqs . 132 and 136 this general frictional 
pressure loss is 

A Pf 2b 2 V M (l-^ 0 ) 

The friction factor X^ is a function of both Reynolds number 
and surface roughness . The value of X^ is obtained from the curve of 
Figure 16, also from Vavra (Ref. 2), for a relative surface roughness 
of 



S- 1 



(lhi) 



where is the peak- to -valley roughness of the surface. Hence, for 
this case 



K 



r 



— j— — 200 b u - in. 
10 



(lte) 



If it is assumed that the rms roughness be one-half K , the surface 
finish must be better than 100 b micro-inches . 

The applicability of this' data to the three dimensional flow 
channels of the energy absorbed is obviously questionable. The losses 
defined below must , at best, be considered as first estimates only. 

The loss coefficient of eq. 129 will be taken as 



C 

u 



2.0 | + 



(1^3) 
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Figure 16 Estimated Channel Wall Friction Factor X £ (Vavra Ref. 2) 



(144) 



and of eq. 130 is 



C . - 2 . 5 § + 



2b X f 



where | is obtained from Figure 15 and X from Figure l6 for 



and 



2b V, 






M3 



u 



-( 



$ UD R 



R 2 *82 b 



R\ 2b s 



2 R 



3 ^3 






2b R V M1 



-( 



R. 



$ a) R, 



2 ^B2^ 

1 



2b 



R 



2 R- W v 



(145) 



(146) 



In eq. 4 the length is given by 



L = 
u 



(L ). + (L ) 

U 1 u o 



where (L ). and (L ) are the lengths of the inner and outer 
v u l u o 

meridional contours, respectively, of the walls of the upper passage 
between stations 2 and 3* Similarly in eq. 133 






<Vi + <Vo 



(147) 



where (L„). and (L„) are the. counterparts of (L.) and (L ) for the 
v Jl'i I'o v l'u v u o 

lower passage between stations 4 and 1. 

5.0 Losses of Angular Momentum in Upper and Lower Passages 

The frictional moments along the upper and lower passages 
were introduced by eqs. l8 and 22 respectively as 



M fu ” X u m E 2 V u2 



(H*8) 



and 



M f* = ^ m R 4 V u4 



(149) 
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No data is presently available for the losses of angular 
momentum in passages of the type needed between the rotor and stator 
of the present absorber. A very crude approximation of these moments 
is obtained, however, by assuming that the shear stresses Y q on the 
walls are merely those acting on a flat plate with a peripheral velocity 
V at R being the velocity outside the boundary layer. This over- 
simplification ignores the effect of the meridional velocity components 
which may cause large increases in the frictional moments . 

For fully turbulent boundary layers along rough walls the 
shear stress is 



Y o = 



c V 
c f 2 u 



(150) 



where c^ is assumed to be invar ient with the Reynolds number of the 
flow. On the element of wall area 2n Rd L of Figure 2 there acts a 
frictional moment 



dM„ = 2n RdL Y R 
i o 

which becomes when using eq. 3 

dM f = rr p c f (R Vj 2 d 1 (151) 

It is assumed that the passages have been designed such that the flow 
is not separated from the walls, and all stresses external to the 
boundary layer are an order of magnitude smaller than stresses along 
the walls, and therefore negligible. This results in the product (RV ) 
being constant outside of the boundary layers and equal to R^V^ , the 
value at the rotor discharge, for the upper passage. Integrating 
eq. 151 from station 2 to station 3 yields 
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(152) 



,(3) 2 

M fu " J (2) " 0 c f < K 2 V “ 

M fu " " 0 C f < E 2 V U 2> 2 [<Vi + <Vo] (153) 



where the inner wall is of length (L^)^ and the outer wall is of 
length (L ) . Hence by eq. 148 



u o 



TT 



\x = 



p c_ R 0 V T(l ). + (L ) 1 
K f 2 u2 L N u i u o.l 



(154) 



m 



Substitution of the mass flow equation and veL ocity relations into 
eq. 154 gives 

1 + § Tan 8 0 f"(L ) . + (L ) 1 
2 L u'i u'oj 



Ni = C f 



2 °R ^2 



(155) 



from eq. 19 






-1-c, 



1 + $ Tan p 2 [(L u ). + (L^J 



2b R ^2 



( 156 ) 



For the lower passage a similar procedure and derivation leads 



to 



= *84 



c f Tan % [ (i A + <Vo] 



2 b 



(157) 



and 






= 1 - 



c f Tan “4 [ ( Vi + ( Vo] 



(158) 



2 ^4 b £ 
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In eqs. 155 through 158, the absolute value of the flow angles 
must be introduced since the loss is a function of absolute deflection. 

For fully turbulent boundary layers, as in the present case, 
the local skin friction coefficients c^ depend on the surface roughness 
ration K /x, where is the peak to average value roughness and x is 
length along the plate. 



c „ “ .004 for K /x = 10" 
f s' 



and 



c “ .006 for K /x = 10 -3 
f s' 

Since there will be flat; irregularities and possible flow -separations , 
it is advisable to use a larger value for c^. A c^ = .016 was chosen 
since it gives best agreement with M-21 test data. 
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AEEENDIX B 



DETAILED ANALYSIS OF THE ENERGY ABSORBER WITH 
THE FLUID CIRCULATING IN THE REVERSE DIRECTION, MODES 3 AND 4 

A separate set of equations must be derived for the reverse flow 
modes. The reason, mentioned earlier, is that the relative flow exit 
angles which are assumed constant are now different. 

This development follows very closely to the previous analysis. 
Therefore, the intermediate steps of this derivation are omitted. 

1. Assumptions 

a. = constant 

a ^ = constant 

The remaining assumptions are identical to the normal flow case. 

2. Analysis 

The geometric relations employed are obtainable from Figures 6 

f 

and 7 * A new flow coefficient § is conveniently defined as 




The torque exerted on the fluid by the rotor when acting in the 
direction of rotation is 

V = " < R 1 V „1 - E 2 V < 2 > 

where the mass flow equations are identical to those for the normal 
case. 
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The frictional moments acting from stations 1 to 4 and stations 



3 to 2 are respectively 



where 



M « = h " K 1 V ul 



M fu = m «3 v u3 



m Ri V | = m R V - M 
4 u4 1 ul fX, 



m R 0 V „ = m R V „ - M _ 

2 u2 3 u3 fu 



The relative flow inlet angles are 



R 
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Tan P 2 = \bx lT 



An 

°3 “ \R ± J K B1 1 / $ 



and 



Tan = 



% (1 + $' Tan p x ) 



' $ 



hi b 
hk b 



The torque is now reduced to 



T r = 2 TT p R. 






2 

(JD T t 



^1 ~ 1 



where 



t h - * + *’ 2 [ Tan Pi - sg s; V T “ “3] 



(3) 



(4) 



(5) 



( 6 ) 



(7) 



(8) 
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Proceeding to the rotor, one can write. 



w 2 = Vr 2 cos ^ 



R 



P 2B " P 2 



Y R = 



P E2 ~ P E1 

£ w 2 

2 1 



and finally 




where 



= ?! " P 2 

, b R 
'k r 2 " R i 



for this case 



n 



A H 



A H 

2 2 
U) R 



In the lower passage 

= p 4Th " p 4 

£ v , 2 

2 m4 



( 10 ) 



( 11 ) 
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and the energy drop is 



II 



1 “ H 4 ^ " 11 w) / R 1\ 2 

2 vr 4 ; 



2 t > 2 

<u R 1 



+ ** pj 

4 



( 12 ) 



+ § 



■2 r°e Ai b R R i\ 2 d - ” m> / E i\ 2 



/_i f gl _R _L' 

12 Vk b 4 b s V ‘ 2 



(£> T “ 2f, i 1 



Nor for the stator 



V. = V. cos i 
4 4s 



— O' — O' 

4B 4 



Y = 



P T4 ~ P T3 
P V 2 
2 v 3 



A 



R, “ R 4 

3 4 



Ap = a_ - or. 
s 3 4 

and the energy drop across the stator is 



H 4- E 4 
U > 2 R 1 2 



. 2 . 

sm i 



,R„, 2 



R H \ 2 



1,2 m (rt) + *' {( sin2i s ) T|2 i« Gr) Tan B i) 



+ $ 



‘2 



• 2 - p 

sm l ,R. 



2 S (r) i w) {(C bp + TaE e il 



tm 2 



(13) 



$ <2 ^R ^BV\ 2 

2cosV 3 S b s V 
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In the upper passage 



C = 



P 2Th " P 2 



u V, 



M2 



and 



H 5- H 2 

cu 2 R 2 
1 



& 2 Wl -n2 S 
2 \R 2 ) i C u VK^J ^ ^ MuHl 



R 



— — Tan or 

S 3 * 



The loss coefficients are now given as 

°u ° 2 -° 5 •+ isr h 

s 



c, = 2.5 S + X r 



h * 



T, 

c « Tan O' - — rr— 

f 3 b B V 



X« = c„ Tan 3 



f ^ b R Si 



The Reynolds numbers ot obtain § and X^ graphically are 



2b , R. 

■W - -V s 0 " 



Sl \ 



*3 S 3 ' 



2b. 



R *» 



S* = “V $ ^ R 



R 1 



The expressions for Y, Y_, and Y remain the same as in the 

K S 



analysis, but the Reynolds numbers are now 



( 14 ) 



} (15) 



(16) 



(17) 



(18) 



previous 
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cosg^ + 



(b R /Ri Tan ^ - Tan $ 2B 



cos £5 



1 



0.425 $ - 1) 



^ 






§’ ou R 1 2 
v 



COSO, + 
2 



b r 

[Tan - Tan 



2 

cos a_ 



,R A \ J4 

°- 425 (^) ( r ^- 1 -° 



(19) 



The equation which is finally solved to obtain $ is 



A L = A H + A H + AH, 

R u s x 

or 

H 1 - H 2 - (H 4 - H 3 ) - (H 3 - H 2 ) - (H 1 - H 4 ) = 0 (20) 

The equation is solved by the program. An iteration technique must be 

1 

carried out until all parameters dependent on $ converge. 
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APPENDIX C 



PROGRAM NORMAL 

Program NORMAL is essentially a modification of the program 
developed by NAEC, 

Program NORMAL solves the equations governing normal circulation, 
Modes 1 and 2, for a given design. The solution is based on an 
iteration technique which is repeated until the flow coefficient 
$ converges to within a specific accuracy. A starting value of § 
is picked as are initial values of the various Reynolds numbers which 
are also $ dependent. The initial iteration step occurs with all 
losses based on standard initial Reynolds numbers. This subsequent 
value of § is then used to improve the estimates of the Reynolds 
numbers and consequential losses. The technique is carried out 
"until $ converges to the proper value. 

If the design parameters and angular speed are input data, the 
program will compute the flow coefficient, torque, velocities, losses, 
energy changes over each component, and the dimensionless torque, 
power, speed, and energy. The flow chart of Table XV depicts this 
process. 

The method of data input is clearly indicated in the program 
itself. The program is felt to be self-explanatory in that comment 
cards have been extensively used. The Fortran symbols employed are 
given in the list of symbols. 

Subroutines have been used liberally. Each subroutine accomplishes 
one individual task such as determining a loss coefficient, or 
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enthalpy change across the stator, etc. The subroutines Fig. 15 and 
Fig, 16 are closest point approximations of the curves depicted in 
Figures 15 and 16 of the analysis. Thses points are given by Block 
data which is called by each subroutine. 
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Table IV C L0^ CHART FOR CO^VVTK R PROGRAMS NORMAL AND REVFRSE 



New 0 and 
Revnolds Mos 
nov; used . 



Data Input: Desiyn c 1ze 

As sien initial values of 
0 and Reynolds numbers 
to bcein iteration. 

N = 0 



N-0 



— jJ 






i 


fi\ vuer 


^ « /V L b' T 




Co- 


r f 






. b, "W 1 


N^/V- | 




R/V'S 


1 

r - V 


i . J 




V 


V ! 


'fl 


•5 l 



i_is I is. 



A/, 



•“b tfyfc 




r l ov Analysis: Calculates and Prints 

Veloc i ties 

Ener^v Differences 

Lo c ses and Loss Co^ r ficient^ 

Non-d i^onM one 1 i zes all pertinent quantities 

Pressure Differences 

Tornue 

Flow Coef fic* ant 
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ANB1 = (ETAML*KB1/KB4*TAN( ALPHA4) )- ( R1R 2**2 ) * ( KB1 /KB2 ) /PH I 
ETA1 = ATAN(TANBl) 

R = BETA1 - BIB 



NA3 = ETAMU*(1. + PH I*TAN ( BETA2 ) ) / ( PHI*KB2/KB3 ) 
PHA3 = AT AN( T ANA3 ) 

= ALP HA 3 - A3B 
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THIS SUBROUTINE ACCOMPLISHES THE FLOW ANALYSIS. ANY MODIFICATIONS TO 
THE DESIRED OUTPUT SHOULD PROBABLY BE MADE IN THIS SUBROUTINE. 

SINCE THE FLOW COEFFICIENT IS ALREADY KNOWN, ANY VALUES IN TERMS OF 
PHI OR GEOMETRIC PARAMETERS CAN NOW BE EVALUATED 
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APPENDIX D 



PROGRAM REVERSE 

Program REVERSE solves the equations for reverse circulation, 
Modes 3 and 4 » The program is very similar to program NORMAL, and 
the flow chart of Table IV is the same. Further discussion would he 
mere repetition of Appendix C. 
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PROGRAM REVERSE 
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CARD SIX: THE RATIO OF THE AVERAGE LENGTH OF THE UPPER PASSAGE TO ROTOR 

BLADE HIEGHTt THE RATIO OF THE AVERAGE LENGHT OF THE LOWER PASSAGE TO STATOR 
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XI = (1.0 - (C0S( IR)**2) )/(COS(BETA2)**2) 

CH1MH2 = 1.0 

BH1MH2 = TAN(BETAl) - ( KB1/KB3) *ETAMU*TAN( ALPHA3 ) 
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